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Abstract 

A simple Almost-Riemannian Structure on a Lie group G is defined by a linear vector 
field (that is an infinitesimal automorphism) and dirn(G) — 1 left-invariant ones. We state 
results about the singular locus, the abnormal extremals and the desingularization of such 
ARS’s, and these results are illustrated by examples on the 2D affine and the Heisenberg 
groups. 

These ARS’s are extended in two ways to homogeneous spaces, and a necessary and suffi¬ 
cient condition for an ARS on a manifold to be equivalent to a general ARS on a homogeneous 
space is stated. 
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1 Introduction 

An almost-Riemannian structure (ARS in short) on an n-dimensional differential manifold can 
be defined, at least locally, by a set of n vector fields, considered as an orthonormal frame, that 
degenerate on some singular set (it can also be defined by a non regular full rank distribution 
and a metric, see Section [2J. This geometry goes back to [20] and (29]. It appears as a part of 
sub-Riemannian geometry, and has aroused some interest, as shown by the recent papers m , 0 , 

HD], HU, m, H31, HU US], M- 

On the other hand the (improperly) so-called linear systems on Lie groups have been studied 
quite thoroughly for ten years, in particular by the authors (see 0 , ®, m1 0 , m and |21j . |22| . 

m, eh mai¬ 
lt is very natural to define an ARS on an n-dimensional Lie group by n left-invariant or affine 
vector fields, the rank of which is equal to n on a proper open and dense subset and that satisfy 
the rank condition. 

As the invariant vector fields can be projected to homogeneous spaces, and affine vector fields 
can be defined on homogeneous spaces, the previous definition of an ARS can be extended to 
homogeneous spaces. In Section [S] it is proved (under some technical assumptions) that an ARS 
that generates a finite dimensional Lie algebra is equivalent (globally or locally, according to 
the technical assumptions) to an ARS on a Lie group or an homogeneous space as previously 
defined (see Theorem a and Corollary [2] of Section [8|. The ARS’s we define on Lie groups and 
homogeneous spaces thus appear as models for an interesting (at least in our opinion) class of 
ARS’s. 
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Excepted in Section [8j devoted to equivalence, we restrict ourselves to what we call simple 
ARS’s, that is ARS’s defined by one linear vector field and n — 1 invariant ones. Notice that the 
famous Grushin plane is of that type, on the Abelian Lie group M 2 . 

After having stated some basic definitions and facts (Section [2| we turn our attention to the 
singular locus, that is the set of points where the vector fields fail to be independant. It is an 
analytic set, but not a subgroup, not even a submanifold, in general. In Section [3] sufficient 
conditions for the singular locus to be a submanifold or a subgroup are stated (Theorems [l] and 
[2j. This locus is very important in what concern the structure of ARS’s, in particular in view 
of a classification. The study is illustrated by many examples, in particular on the Heisenberg 
group, which show that the ARS’s on Lie groups, even the simple ones, go much farther than 
the 2D or 3D generic case (see 0, C3|). 

The Hamiltonian equations of the PMP are computed in Section [4] These computations are 
rather standard, but more complicated than in the invariant case, because the Hamiltonians can 
here depend on the state. They allow to get a complete characterization of the abnormals, it is 
the object of Theorem [3] in Section [5] 

Before providing examples we show that simple ARS’s can be globally desingularized in a 
standard way: we establish in Section [6] a relation between a simple ARS and a codimension one 
regular sub-Riemannian structure on the semi-direct product of the Lie group with the real line. 

In Section [7] two examples are computed, the first one on the two-dimensional solvable group, 
and the other one on the Heisenberg group. 

The computation of some diffentials and power series is postponed to the appendix. 

In order to avoid to lengthen the paper, we have chosen to not recall the basic facts about 
optimal control and sub-Riemannian geometry. The reader is referred to textbooks as [26j, [I], 

0 , m- 

2 Basic definitions 

2.1 Linear vector fields 

In this section the definition of linear vector fields and some of their properties are recalled. More 
details can found in [22]. 

Let G be a connected Lie group and g its Lie algebra (the set of left-invariant vector fields, 
identified with the tangent space at the identity). A vector field on G is said to be linear if its 
flow is a one parameter group of automorphisms. Notice that a linear vector field is consequently 
analytic and complete. 

The flow of a linear vector field X will be denoted by (<£g)tgR. 

The following characterizations will be useful in the sequel. 

Characterization of linear vector fields 

Let X be a vector field on a connected Lie group G. The following conditions are equivalent: 

1. X is linear; 

2. X belongs to the normalizer of g in the algebra V U (G) of analytic vector fields of G, that is 

VTeg [X,Y]eg (1) 


and verifies X(e) = 0; 


Vg ,g G G Xggi — TL g .X g i +TR g i.X g 


3. X verifies 


(2) 


According to ([!]) one can associate to a given linear vector field X the derivation D of g 
defined by: 

VY £ g DY = -[X,Y], 

that is D = —ad(A’). The minus sign in this definition comes from the formula [Ax, b] = — Ab in 
MY It also enables to avoid a minus sign in the useful formula: 

VY eg, Vi e M <£>4 (exp Y) = exp (e tD Y). (3) 

In the case where this derivation is inner, that is D = — ad(X) for some left-invariant vector 
field X on G, the linear field splits into X = X + X* X, where X stands for the diffeomorphism 
g e G i— > 1(g) = g- 1 ■ Thus X is the sum of the left-invariant vector field X and the right- 
invariant one X*X. The flow of such a vector field is given by 

Vt(g) = exp(-tX)gexp(tX) (4) 

An affine vector field is a element of the normalize!' 94 of g in V U> (G), that is 

94 = norm v * (G) g = {F £ Y"(G); VY e g, [ F , Y] € g}, 

so that an affine vector field is linear if and only if it vanishes at the identity. 

It can be shown (see [9] or [22]) that an affine vector field can be uniquely decomposed into 
a sum 

F = X + Z 

where X is linear and Z right-invariant. 

2.2 Almost-Riemannian structures 

For all that concern general sub-Riemannian geometry, including almost-Riemannian one, the 
reader is referred to [Tj. 

Definition 1 An almost-Riemannian structure on a smooth n-dimensional manifold M is a 
triple (E, /, (.,.}) where: 

1. E is a rank n vector bundle on M; 

2. f : E i —> TM is a morphism of vector bundles; 

3. (E, (.,.)) is an Euclidean bundle, that is (.,.) is an inner product on the fiber E q of E, 
smoothly varying w.r.t. q; 

assumed to satisfy the following properties: 

(i) The set of points q £ M such that the restriction of f to E q is onto is a proper open and 
dense subset of M; 

(ii) The modulus S of vector fields of M, defined as the image by f of the modulus of smooth 
sections of E satisfies the rank condition. 

Remarks 

1. The set of points of M where the rank of f(E q ) = S q is less than n is called the singular 
locus of the ARS and will be denoted by Z in the sequel. If M \ Z was not required to be 
proper (see (i)), that is if Z could be empty, then the structure could be Riemannian. 


2. This structure is trivializable if (E, (.,.)) is isomorphic to the trivial Euclidean bundle 
M x M n . In that case we can choose an orthonormal basis (ei,..., e n ) on M n , and define n 
vector fields on M by fi(q) = f(q , e*), i = 1,..., n. The set (/i, ..., f n ) is an orthonormal 
frame on M \ Z. 

Norm 

The almost-Riemannian norm on S q is defined by 

|M| = min{||u|| ; u G E q and f{u) = v}. 

Notice that thanks to the linearity of / on the fibers, the minimum is well defined. In the 
trivializable case one has: 


u = mm 


N 


n 

u l/lO?) H-1- U n f n (q) = V 

1 


Admissible curves 

A curve 7 : [0, T] 1 —> M is admissible if it is Lipschitz continuous (w.r.t the differential 

structure of M) and if there exists a measurable bounded function 1 1 —> u(t ) from [0,T] into E 
such that 7 (t) = f(u(t)) for a.e. t €E [0,T], Its length is: 

Kl )=[ IItWII dt 

Jo 

In the trivializable case, assuming that ||w(f)|| realizes the minimum of || 7 (t)||, we get 


2 ( 7 ) 



The almost-Riemannian distance is define as the inhmum of the lengths of the admissible curves, 
and it is a well-known fact that this distance is finite, continuous, and induces the manifold 
topology (the rank condition is here essential). 


2.3 Simple ARS’s 

In that paper we mainly deal with what we call simple ARS’s. They are ARS’s defined on a 
connected Lie group G by a set of n vector fields 

{A,yi,...,y n -!} 


where 


(i) X is linear; 

(iii) Yi ,..., Y n ~ 1 are left-invariant; 

(iii) n = dim G and the rank of X, Y\ __, Y n -\ is full on a non empty subset of G; 

(iv) the set { X , Yi,..., Y n _ 1 } satisfies the rank condition. 

The set where the rank of X, Y[.... ,Y n -1 is not full will be refered to as the singular locus 
and denoted by Z. 


Remarks. 



1. The singular locus Z is an analytic subset of G. By Assumption (iii) it is not equal to 
G, and by analycity its interior is empty. On the other hand X(e) =0 and it cannot be 
empty. Finally G\Z is an open, dense and proper subset of G. 

2. The rank condition implies 


CA{X, Yi,..., Yn-x} = M-T 0 g. (5) 

Necessary conditions for the rank condition 

Let us denote by A the vector subspace of 0 generated by Y \,..., Y n —\, and by D, as usual, 
the derivation asociated to X. If [A, A] C A and D( A) C A, then the Lie algebra generated 
by X ,Yi,... ,Y n _ i is equal to MA © A. But the rank of that Lie algebra is not full on Z. 
Consequently the rank condition implies that at least one of the following conditions hold: 

(i) [A, A] g A 

(ii) D( A) g A 

Notice that in all cases, the full rank is obtained after one step. 

Example. The Grushin plane 

The Grushin’s metric on the plane corresponds to the frame X = (0,xi) T and Y\ = (1,0) T , 
that is to the control system: 

x = vX + U\Y\ or, in coordinates, < Xl U1 

\ X 2 = vxi 

The state space M 2 is an Abelian Lie group, the field X is linear and the held Y\ invariant, 
so that the Grushin plane is a simple ARS. The singular locus is here the line {aq = 0}. 

2.4 Notations 

In the sequel the following notations will be used: 

1. Let Y 6 g. Then Y g will stand for TL g .Y, where L g is the left translation by g, and TL g 
its differential. 

2. To the linear vector held X we associate F(g) = TL g -i.X g € 0 . We also write F g for F(g) 
when it is more convenient. 


3 Singular locus 

3.1 Introduction 


This section is devoted to the properties of the singular locus Z of the simple ARS determined 
by the orthonormal frame {X,Y\,... ,Y n -\}. The metric does not matter here and instead of 
Y \,..., y n _i we will mainly consider the subspace A of 0 they generate. 

In the Abelian simply connected case, that is when G = M n , the singular locus is easily seen 
to be a codimension one subspace of M n . It is natural to ask whether Z is a subgroup of G in 
general, but the examples of Section 3.5, in particular on the group Heisenberg, show that it is 
not the case. Actually Z is not even a manifold in general, but only an analytic set. 

By definition 


Z = {g eG / rank (X (g), Y 1 (g ), Y 2 (g ),.. • , 1^-1 (g)) =n - 1}. 



However we need to characterize Z in a more handable way. The hyperplane A can be defined 
as the kernel of a one-form oj to which we can associate a left-invariant one-form likewise denoted 
by u. 

We recall from the end of the previous section that the mapping F from G to g is defined by 

F(g) = TL g -i .X g 

The power series expansion of F and its differential are computed in the appendix (Propositions 
[4} [ 5 ] and Corollary [ 3 ]). Thanks to ui and F, the singular locus can be defined in a form that allows 
differential calculus: 


g & Z X g G A g <*=/- F g £ A <^=^> (w, F g ) — 0. 

In the sequel we denote by if the analytic function from G to M defined by ip(g) = (oj, Fg). The 
singular locus is consequently the set of zeros of i/j. 

3.2 A sufficient condition for Z to be a submanifold 

In this section it is shown that Z is a submanifold (embedded and analytic) as soon as A is a 
subalgebra of g. 

Lemma 1 Let us assume that A is a subalgebra of Q. Then g can be decomposed as follows: 

0 = MTn © IT © M£ (6) 


where 

1. IT is a (n — 2)-dimensional subalgebra; 

2. A = IT © M£; 

3. D~ 1 (A) = My n © IT; 

4- (u,Y n ) = 1 and D£ = Y n (mod A). 

Proof As A is a subalgebra, the rank condition implies that D*u = uj o D is different from uj 
and does not vanish (otherwise we would have DA C A and the rank condition would not hold, 
see Section [ 2 ]). Consequently D~ l A = ker (D*uj) is a (n — 1 (-dimensional subspace of g different 
from A, and 

IT = A n D -1 A 

is a (n— 2)-dimensional subspace of g which turns out to be a subalgebra. Indeed let X\. X 2 G IT. 
Since A is a subalgebra we have [Aj, Xf € A. By definition of IT, DX\ and DX 2 belong also 
to A so that 

D[X 1 ,X 2 ] = [DX u X 2 ] + [X]. dx 2 ] G a 

In conclusion [Xi, X 2 ] belongs to D~ 1 A hence to IT which is thus a subalgebra of g. 

Let Y n £ g such that D -1 (A) = MT re © IT. As Y n (f A, ( ui,Y n ) does not vanish and can be 
assumed to be equal to 1. 

To finish DA A hence Df A, and Df = aY n (mod A) with a / 0. Up to a modification 
of £ we can choose a = 1. 

□ 

Theorem 1 If A is a subalgebra of g then the singular locus Z is an analytic, embedded, codi¬ 
mension one submanifold of G. 

Its tangent space at the identity is ker (D*u) = D~ 1 A. 


Proof. 

As if is analytic, we have only to show that its differential does not vanish on Z. According 
to Corollary [3] of the appendix it is given by: 

V<7 G G, YY G o T g if.Y g = (u, DY + [F g , Y]) 

Consider £ as defined in Lemma [l] Since £ G A, F g G A if g G Z, and A is a subalgebra, we 
have [Fg, (]gA for all g G Z. On the other hand D £ ^ A and: 

Vff G A Tglf.fg = (u, Df + [Fg, £] ) = 1. 

This shows that if is a submersion at all points of Z which is thus a codimension one submanifold 
of G. Moreover F e = 0 so that 

VY G g T e if.Y = (cj,DY) = (D*lo,Y). 


□ 

Remark. The previous formula holds at every point where F vanishes, that is on the set of 
singularities of X . 

3.3 Sufficient conditions for Z to be a subgroup 

In order to state sufficient conditions for Z to be a subgroup we begin by a lemma which is not 
easily checkable in practice but which is essential in the proofs of the forthcoming Proposition [l] 
and Theorem [2] 

Lemma 2 The singular locus Z is a subgroup of G if and only if 

Vg, g G Z (lo, Ad{g')F g ) = 0 (7) 

that is if and only if for all g,g' G Z, Ad(g')F g G A. 

Proof. Let g,g' G G. Then, according to Formula ([2]) of Section [2] 

9 ) = TL^g'g^xXg'g 

= TLg-lTL(gl}-l(TLglXg + TRgXgl ) 

= F(g)+Ad(g- 1 )F(g / ). 

Let us assume that g,g' G Z, hence that F(g) and F(g') belong to A. Then 

g g G Z F(g'g) gA^ Ad (s" 1 W) G A. 

In the same way we get 0 = Flgg -1 ) = F^g -1 ) + Ad (g)F(g), so that F(g -1 ) = —Ad {g)F g , and 
g^ 1 belongs to Z if and only if Ad {g)Fy G A. 

Finally we see that Z is a subgroup of G if and only if Formula ([Tj) holds. 

□ 

Proposition 1 The following conditions are equivalent: (i) A is an ideal of g, (ii) the derived 
algebra T >1 g is included in A, (in) the left-invariant one-form oj is closed. 

When these conditions are satisfied the singidar locus Z is a Lie subgroup of G whose Lie 
algebra 3 is equal to ker(Z?*w) = D^ 1 A. 


Proof. 

The first thing to notice is that Z being closed is a Lie subgroup of G as soon as it is a 
subgroup. 

It is well-known that for a left-invariant one-form u the Maurer-Cartan formula writes 

MX, Y E g du(X, Y ) = -u([X, T]) 

so that du vanishes at e (hence everywhere since it is left-invariant) if and only if the derived 
subalgebra T> l g is included in ker(cu). If Z E g is such that g = A +MZ then: 

V 1 g = [A,A] + [A,Z\. 

It is clear that V 1 g C A if A is an ideal of g. Conversely if A is not an ideal then either A is 
not a subalgebra or [A, Z] is not included in A. In both cases T> l g is not included in A. 

Let us assume A to be an ideal of g, and let g , g' E Z. The group G being connected, there 
exist X±, ..., Xk E g such that g = exp(Xi)... exp(Xfc), and: 

Ad (g) = Ad(exp(Xi)... exp(X fc )) 

= Ad(exp(Ai))... Ad(exp(X fc )) 

_ e ad(Xi) _ _ _ e ad(x k ) 

Each of the e a( -l(X) : s senc [ s into A. Indeed for any Leg 

+OO 1 

e ad(A)y = y + V" l a d j (XPY = Y (mod V 1 q) 

< n I 

3 =1 J ' 

and e a d(^i)y g A if and only if Y E A. 

This shows that Ad(cjr) sends A into A, and since F(g') E A we obtain Ad (g)F(g') E A. 
According to Lemma [2] the singular locus Z is a subgroup of G. 

The Lie algebra of Z is its tangent space at e. Since A is an ideal we know by Theorem [l] 
that 3 = T e Z = ker(D*u). 

□ 

Theorem 2 If the Lie algebra g is solvable, and A is a subalgebra of g, then the singular locus 
Z is a codimension one subgroup of G whose Lie algebra is 3 = D~ 1 A. 

Proof. 

Notice first that according to Theorem [l] the singular locus Z is a codimension one subman¬ 
ifold of G whose tangent space at the identity is 3 = D~ 1 A. Therefore we have only to prove 
that Z is a subgroup. 

1. To begin with let 6 be a codimension one subalgebra of a nilpotent Lie algebra I), and let 
us show that 5 is an ideal of I). Let Z E I) such that fi = (5 + MZ; then its derived algebra is 
V 1 fi = [5, <5] + [5, Z], If 5 is a subalgebra but not an ideal of fi then [d, <5] C 6 but P 1 !) (2 5, 
so that \5,Z\ is not included in S. Therefore there exists X E <5 such that ad (X)Z = aZ 
(mod 5), for some real number 0 / 0. We can choose a = 1 so that 

VA: > 1 ad k (X)Z = Z (mod 5) 


which proves that f) is not nilpotent, a contradiction. 


2 . Let now g and A as in the statement. If A is an ideal of g then the result comes from 

Proposition [I] If not D 1 g is not included in A and there exists Z £ P 1 g such that 

g = MZ 0 A. Moreover the rank condition implies that D( A) (Z A. 

3 . Let r) be the largest nilpotent ideal of g (see D 3 ], § 4 , n° 4 ). It is known that for any deriva¬ 
tion d of the solvable algebra g, the inclusion d(g) C 77 holds, so that 77 is a characteristic 
ideal of g (see HZ], § 5 , Proposition 6 and its corollary). Moreover P 1 g being a nilpotent 
ideal of g is included in r/, hence Z £ 77 and g = 77 + A. Consequently 

dim(r/ n A) = dim(?y) + dim(A) — dim(g) = dim(77) — 1 , 

which proves that 77 n A is a codimension 1 subalgebra of 77. Since 77 is nilpotent, and 
according to Item 1 , 77 n A is an ideal of 77. Actually 77 n A turns out to be an ideal of g. 
Indeed any Y £ g writes Y = X + aZ where A £ A and a £ M. Since Z £ 77 and 77 n A is 
an ideal of 77 it is clear that ad(Z) sends 77 flA into 77 n A. The same is true for ad(X) 
because X belongs to the algebra A and 77 is an ideal of g. Finally 77 n A is invariant for 
all inner derivations, hence an ideal of g. 

4 . Let us show that D~ 1 A is a subalgebra. Firstly we have D(D~ 1 A) C D(g) C 77 and 

D(D~ 1 A) C A by definition, so that D(D~ 1 A) C 77 n A. Then for all X\, X 2 in D~ 1 A 

we get DXi £ 77 n A, hence [DXi,Xj\ £ 77 n A for i, j = 1 , 2 . Consequently 

D[X 1 ,X 2 ] = [dx u x 2 ] + [x u dx 2 ] £ 7? n a 

and [X\, X 2 ] £ D~ 1 A, which finishes the proof. 

5 . In this item we show that F(g) £ 77 for all g £ G. According to Proposition [5] of the 
appendix 

MY £ g, F(exptF) = y(-l) fc - 1 7 -ad (fe - 1) (y)T>y £ ?? 

re! 

k> 1 

because DY belongs to the ideal 77. Secondly 

F(g exp tY) = F(exptF) + e~ tad ^ Y) F{g) 

belongs to 77 as soon as F(exp tY) and F(g) both belong to 77. Since g = exp(Ti)... exp (Yf.) 
for some Yf,..., Y^ £ g the result is obtained by induction. 

6. To finish the proof of the theorem let g, g' in Z. On the one hand F(g) belongs to the 
ideal 77 n A. On the other one g' = exp(Ai)... exp(Afc) for some X\, ..., Xj~ £ g, hence 

Ad ( ff ') = e ad(Xl) ...e ad ( Xfc ) 

and Ad (g')F(g) belongs to 77 0 A. The singular locus Z is a subgroup of G according to 
Lemma | 2 ] 


Remarks. 


□ 


1. The beginning of the proof shows that when g is nilpotent then A is an ideal as soon as it 
is a subalgebra. In the solvable case it is easy to exhibit examples where A is a subalgebra 
but not an ideal (see Example 3.5.2 on Aff + (2)). 



2. The same is true on semi-simple Lie groups since a codimension one subalgebra cannot be 
an ideal. But apart from the fact that codimension one subalgebras are not very common 
in semi-simple Lie algebras, that condition is not here sufficient for Z to be a group (see 


Example 3.5.4 on SL2). 


There is another case where we can assert that Z is a subgroup, it is when it is equal to the 
set of singularities of X. 

Proposition 2 If Z is equal to the set Zx of singularities of X then it is a closed Lie subgroup 
of G whose Lie algebra is equal to ker(D). 

Proof. The set Zx is obviously a closed subgroup of G, since 

Z x = f}{geG; <p t (g) = g}, 

te K 

where (tpt)te .R stands for the flow of X. Moreover 
Vs E M exp(sT) E Zx 

<*=>■ Vs G M, Vt E M, exp(sT) = ^(exp(sT)) = exp (e tD sY) 

<=^ Vt E R, Y = e w Y 
^ DY = 0, 


which shows that the Lie algebra of Zx is kerZL 


□ 


This case is illustrated in Examples 3.5.2 and 3.5.3 


3.4 Necessary conditions for Z to be a subgroup 

When the one form D*u does not vanish, the singular locus is locally, around the identity, a 
submanifold, whose tangent space at e is T e Z = ker(D*cu) = D~ 1 A. One could expect Z to be 
a subgroup if D~ l A is a subalgebra, but this statement is wrong, even locally, and even if ZW 1 A 
is an ideal of g (see the examples of Section 3.5 in particular Example 3.5.3.1). 

In order to state a necessary condition for Z to be a subgroup, and a locally sufficient one, 
we will use the following stronger algebraic condition: 

Let f} be a subalgebra of g. It will be said to satisfy the condition ( HZ ) if 

V Y E f), V m > 0, V Z\,... , Z m G f), (w, ad(Zi)... ad (Z m )DY) = 0. 


Proposition 3 If Z is a group, then it is a Lie subgroup of G whose Lie algebra 3 satisfies 
Condition (HZ). 


Proof. The set Z being closed is a Lie subgroup as soon as it is a group. 

For all Zi,... , Z m G 3 and for all t \,..., t m GR the point exp(tiZi)... exp (t m Z m ) belongs 
to Z. Since Ad(exp(iiZi)... exp (t m Z m )) = e had(Zi) _ _ _ e t m ad(z m )^ an( j according to Lemma [2J 
we get for all g G Z 

(u, e L ad (^i)... e tmad(z m ) F ^ = a 

Derivating this equality at iq = 0,..., t m = 0, and thanks to the linearity of ad (Zi), we obtain 


(cu, ad(Zi)... &d(Z m )F(g)} = 0. 






To finish let g = exp(tT), where Y £ 3. According to Formula ( |14[ ) (Proposition [d| Appendix), 
we have 

4- F(exptY) = DY, 
di|t=o 

which gives the desired equality: 


(w, ad(Zi)... ad (Z m )DY) = 0. 


□ 


Remarks 

1. Proposition [3] is a local result: actually it holds as soon as Z is locally, in a neighbourhood 
of the identity, equal to some Lie subgroup of G. 

2. The successive derivatives of F(ex.p(tY)) do not provide more information since some of 
the Zi s can be chosen equal to Y. 

3. Condition (HZ) is trivially satisfied only if A is an ideal. 

In order to investigate the converse to Proposition [3] we will consider a subalgebra h of 0, 
assumed to satisfy Condition (HZ), and we will associate to f) the following subspace of A: 

5h = {X £ A; V m > 0, V Z \,..., Z m £ f) ad(Zi)... ad (Z m )X £ A}. 

Notice that Sh depends on f) and, thanks to Condition (HZ), that D f) C Sh- It is also clear that 
5h is ad(Z)-invariant for all Z £ f). 

Lemma 3 Let f) be a subalgebra of q that satisfies Condition (HZ), and let H be the connected 
subgroup generated by f). Then 

1. The Lie subgroup H is included in Z. 

2. Let g £ Z. If F(g) £ Sh, then gH is included in Z. 

Proof. 

1. Let Y £ fj and k > 1. Condition (HZ) with m. = k — 1 and Zj = Y for i = 1,..., k — 1 is: 

ad ( k ~ 1 \Y)DY £ S h C A 

Consequently for all t £ M 

jk 

F(exptT) = ^2(~l) k - 1 -ad^- 1 \Y)DY £ 6 h C A 
fc>i 

and exp (tY) belongs to Z for all t £ M. 

2. Now let g £ Z such that F(g) £ 5} l . For any Y £ f) we have 

(w, F(g exptY)} = (u, F(exptY) + ^F(g)) 

= {u, F(exp tY)) + Efc>0 T ad(fc) ( Y ) F (a)) 

= 0 

because on the one hand the first term vanishes according to the first item, and on the other 
one ad ( k \Y)F(g) belongs to A since F(g) belongs to Sh- Notice moreover that F(gexptY) 
belongs to Sh, since Sh is ad(Z)-invariant for all Z £ f). 



3. Let Z i,..., Z m in f) and t\,..., t m in M. According to Item 2 we obtain by induction that 
exp (tf,Zk)... exp(tiZi) belongs to Z for k = 1,..., m, but moreover that 
F(exp(tkZk) ■ ■. exp(iiZi)) belongs to 5h, allowing the induction. 

This shows that H is included in Z. 

4. To finish let g E Z such that F(g) E 5h- We get gH C Z by the same reasoning as in Item 
3 . 


□ 

Corollary 1 If D*lo is not zero and if 3 = ker [D*oj) is a subalgebra of g that satisfies Condition 
(HZ), then the connected subgroup Gr(]f) generated by 3 is included in Z. 

Moreover there exists a neighbourhood V of the identity such that: 

Z n V = Gr(i ) n V. 


Proof 

The first part is immediate from Lemma [3] For the second one the condition D*lo 7 ^ 0 
implies that Z is locally, in a neighbourhood V of e, a (n — l)-dimensional submanifold (recall 
that Z = {?/> = 0} and T e fj = D*iv). This submanifold containing the (n — l)-dimensional Lie 
group Gr( 3 ), is equal to that last within V. 

□ 


The examples 3.5.3.1 on the group Heisenberg (see Section 3.5.3) show that the assumptions 
of Corollary [l] are not sufficient to know what happens far from the identity: they imply neither 
that is regular in Z nor that one has F(g) E 5* for all the points g E Z, even in the connected 
component of Z. 


3.5 Examples of singular loci 
3.5.1 Abelian groups 

Consider the ARS on M n defined by {A, b±, ..., 6 n _i) where A is a real matrix of order n and 
b \,..., b n -i are n — 1 linearly independent constant vector fields . It is straightforward to show 
that the singular locus is a codimension one subspace of M n . 


3.5.2 The 2D affine group 

Let G be the connected component of e in the 2-dimensional affine group: 


G = AfU( 2) = {(o 1 ); ( x i v) £ x M| . 


Its Lie algebra g = off(2) is solvable, generated by 


X = 


1 0 
0 0 


and Y = 


0 1 
0 0 


with [ X , Y ] = XY — YX = Y, and identified with the set of left-invariant vector fields which is 
therefore generated by 


gX = 


x 0 
0 0 


and gY = 


0 x 
0 0 


where g = 


x 

0 


y 



is equal 


All the derivations are inner and the one whose matrix in the basis (A, Y) is 


0 0 
a b 


to D = —ad (aY — bX). The associated linear vector held X is given at the point g by 




0 a(x — 1 ) + by 

0 0 


Let us consider the simple ARS defined by X and Y\ = aX + j3Y = 


ax fix 

0 0 


In order that 


X and Y\ be independant on some subset the constant a must not vanish. 

Then DY\ = aaY + b[3Y and the rank condition is equivalent to era + /3b ^ 0. 

As x > 0 the first component of Y\. that is ax , never vanishes, so that Y\ and X are colinear 
if and only if X vanishes. Finally we get 

Z = {g e G\ X g = 0} = {(x, y ) £ X M; a(x — 1) + by = 0}. 

According to Proposition [2] it is a closed Lie subgroup of G. Another proof comes from the fact 
that A is always a subalgebra since it is one dimensional. As the group Aff + ( 2) is solvable we 
know by Theorem [2] that Z should be a subgroup. 


3.5.3 The Heisenberg group 

Let G be the simply connected Heisenberg Lie group of dimension three 


G = 


r a 

X 

A 

0 

1 

y 

l \o 

0 

i J 


Its Lie Algebra 0 is generated by 



where [X, Y] = XY — YX = Z and the other brackets vanish. As left-invariant vector fields 
they write in natural coordinates: 


X = — 
dx ’ 


v d d 

~ dy +X dz' 


z =k- 


The derivations of g are the endomorphisms D whose matrix in the basis (X, Y, Z ) has the form: 


fab 0 \ 

D = I c d 0 J 
\e / a + dj 

and the associated linear vector held is: 

d d 1 1 (9 

X(g) = (ax + 6y)— + (cx + dy)— + (ex + fy + (a + d)z + -cx 2 + T^V 2 )-^- 

More details can be found in [22] . 

Let us now consider a simple ARS on the group Heisenberg, defined by a 2-dimensional 
subspace A of g and a derivation D. We get two very different behaviour depending on whether 
A is a subalgebra or not. 

One can also notice that the number of connected components of G\Z ranges from 1 to 4. 


1. A is a subalgebra. According to Theorem [2] the singular locus Z is a subgroup of G. Up 
to an automorphism of g we can assume that A = Span{X, Z}. Then the rank condition 
is satisfied if and only if c 7 ^ 0 (in the matrix of D). A particular case is when Z is exactly 
the set of singularities of X. One example is obtained by 

/0 0 0 \ 

D = 1 0 0 

\e 0 0 / 

but is far from being the only one. Here X = x-j^ + (ex + ^x 2 ) -jp and it is easily seen that 
Z = {X = 0} = {x = 0}. 

2. A is not a subalgebra. We can assume without lost of generality that A = Span{A, Y}. 
It is easy to see that the ARS is well defined for any derivation different from zero. Indeed 
the rank condition is satisfied and 

Z = {ex + fy + (a + d)z - ^ cx 2 + ^ by 2 - dxy = 0} 

is equal to M 3 if and only if all the coefficients of D vanish. 

It is also clear that the loci defined by these quadratic forms need not be subgroups, not 
even submanifolds. Let us exhibit some particular examples. 

Example 3.5.3.1 Let D have the following form: 

(a b 0 \ 

D = I c —a 0 I . 

\0 1 0 / 

In the {A, Y, Z} basis we have ui = (0,0,1) and D*lo = (0,1,0). 

The one form D*uj does not vanish, and its kernel, which is also the tangent space to Z at 
the origin, is the subalgebra 3 = D~ 1 A = Span{X, Z}. However this fact is not sufficient 
for Z to be a subgroup, even locally. 

If c 7 ^ 0 the condition (HZ) is not satisfied, and Z = {y — \cx 2 + \by 2 + axy = 0} is 
not a local subgroup around the origin: it is straightforward to check that (x,y, z)~ 1 = 
(—x, —y, —z + xy) does not belong to Z in general when (x, y, z) E Z. 

If c = 0 then 3 satisfies Condition (HZ) so that we know by Corollary [I] that Z contains 
the group generated by 3 , that is Gr( 3 ) = {y = 0}. However 

^ = {y( 1 + ^ by + ax) = 0} 

and Z reduces to {y = 0} if and only if a = b = 0. Otherwise Z is the union of the plane 
{y = 0 } with another plane which may intersect it (if a / 0 ) or be parallel to it (if a = 0 ). 

Let us consider the case a = 1. 

At the points g = (—1,0, z), that belong to Gr( 3 ), the function ip is singular, that is T g ip 
vanishes. 

At the points g = (—1 — \by, y, z) with y ^ 0 , it is clear that gGr( 3 ), the translation by g 
of the connected group generated by 3 , is not included in Z. 


Example 3.5.3.2. Degenerated case 


Consider 

/0 b 0 \ 

D = I c 0 0 I with 6 > 0 and c < 0 

\0 0 0 / 

Here D _1 A = g and the condition ( HZ ) is not satisfied. 

If b > 0 and c < 0, then Z is the codimension 2 subgroup {x = y = 0} of g. 

If b and c are both positive (or negative), then Z = {bx = Aby} is the union of two secant 
planes. 

If b ^ 0 and c = 0, then Z is the codimension 1 subgroup {y = 0} of g. 


Example 3.5.3.3. Tangential case 
Let 

/a b 0 \ 

D = I c d 0 I with a + d ^ 0. 

\0 0 a + d) 

The one forms uj = (() 0 l) and D*u = (0 0 a + d) define the same subspace A = 
D~ 1 A of g. As the differential at the origin of the defining function of Z. that is T e T = D*u 
does not vanish, the singular locus is a submanifold in a neighbourhood of e, and e is a 
tangential point, that is A = T e Z. 

Unlike the generic case (see [13]) the tangential points need not be isolated. Consider 


D = 



The singular locus is Z = {z = x 2 + xy}, and A is tangent to Z along the parabola 
{z = — x 2 } contained in the plane {y = — 2x}. 


3.5.4 The special linear group SL( 2;M) 


Example 3.5.4.1. Let G = <SX(2;M) be the order 2 special linear group. Its Lie algebra 
g = sl(2;M) is the set of matrices with trace zero. The usual basis: 


of sl( 2 ; M) verifies 



[H,X] 


0 ), 

X = 

0 !\ 

-v’ 

°y 

= 2 X, 

[H,Y] 

= - 2 T, 



[X, Y] = H. 


Consider the simple ARS defined by A = Span{Lf, X} and the derivation D = — ad(T) (sl(2;M) 
being semi-simple, all the derivations are inner). We get D~ 1 A = Span{A, Y} which is not a 
subalgebra. Consequently the singular locus cannot be a codimension one subgroup of SL( 2; M) 
despite the fact that A is a subalgebra of sl(2;M). This shows that Theorem |2] does not apply 
to general groups. 


Let us compute the singular locus. At the point g 


a b 
c d 


with ad — be 


1 , the linear 


b 0 
d — a —b 


A straightforward 


vector field associated to D = — ad(y) is X g = gY — Yg = 

computation shows that the vectors X g . H (j , X ;) are linearly dependent if and only if a = ±1, in 
other words that 


Z = {a = ± 1}. 


Thus the singular locus Z is a submanifold of SL( 2;M) but nor a subgroup neither a connected 
set. 

Notice that the ARS is well defined: indeed Z ^ G and the rank condition is satisfied. 

Example 3.5.4.2. Another example can be obtained by switching the roles of A and T>~ 1 A. 
More accurately let A = Span{A, Y}, it is not a subalgebra, and D = —ad(A). We get easily 
D~ 1 A = SpanjiL, A'}. It is a subalgebra of s((2;M) that satisfies Condition (HZ), hence the 
singular locus contains the connected subgroup generated by D~ 1 A. However this locus turns 
out to be Z = {cd = 0}, which has three connected components: 



Ci is a connected subgroup, Ci U C is a non connected subgroup, but Z = C\ U C2 U C3 is 
clearly not a subgroup of SL( 2;M). 

4 Hamiltonian equations 

The aim of this section is to state the Hamiltonian equations of the PMP applied to ARS’s on 
Lie groups. We follow the same lines as in the invariant case (see for instance 0 ) but here the 
Hamiltonian depends on the point g £ G which entails some more complications (see also | 28 j). 

4.1 The canonical symplectic structure of T*G 

All the material of this subsection is standard and can be found in the previous references. 

As usual the cotangent bundle T*G is identified to g* x G by 

4- : 0 * x G —> T*G 

(A ,g) 1 —* A s = Ao TL g -i G T*G, 

so that the equality ( \ g ,Y g ) = (A ,Y) holds for all Y 6 g (recall that Y g stands for TL g .Y). 

The projection from T*G to G is denoted by n in what follows. 

The tautological one-form s on T*G being defined by (s^, () = (A, n*£), its pullback by <t> is 
the one-form 4>*s defined on g* x G by (here X 6 g): 


{&s^ g) ,(ti,X g )) = ( X,X). 


In this setting, the symplectic form d&*s = &*ds is characterized by 

**ds {a>9) ((e, X g ), (77, Y g )) = <e, Y) - (rj, X) - (A, [X, Y}) . 

To finish this review let h G C°°(T*G ) and let h be the Hamiltonian vector held associated 
to h by dh = —i^ds. If H is equal to h o $ and H is defined by h = it is easily verified 

that H is the Hamiltonian vector held associated to H, i.e. 

dH = -i g $*ds. 


4.2 Computation of the Hamiltonian vector fields 

In what follows we set a = &*ds. 

Let H be an Hamiltonian on g* X G and H the associated Hamiltonian vector held. We can 
identify the tangent space to g* x G at the point (A, g) with g* x T g G, and write H = (£,X g ) at 
this point. Then for all ( 17, Y g ) E g* X T g G we get on the one hand 


dH( n ,Y g ) 


/dH 

\-d\’ V 


+ / d JL Y 

+ w 3 


and on the other hand 


dH( V ,Y g ) =-a {X!g) ((Z,X g )(r ] ,Y g )) 

= -{£,Y) + {ri,X) + { A, [X,Y]). 


1. Setting Y = 0 we get X 

2. Setting 77 = 0 we get 


^-(A,s) e (g*)* = g. 

(^ , ^) = -(C^ > + (A,[x,y ]> 


hence 


so that 


Summarizing we obtain 


(t,Y) 


= ((ad(X))*\,Y)-((TL g r^L 

£ = (ad(X))*\ -(TL g )* — . 



where 



- (dH\ 8 (( ,,dH. 

E ={ax),^ + {{^ 

dH 

should be understood as TL g ——. 

o A 


A - (TL g )* 


dH A 
dg ) 


d_ 

ax' 


4.3 Left-invariant and linear vector fields 

The Hamiltonian associated to a left-invariant vector held Y. that is H = (A, Y), does not depend 
on g, so that the corresponding Hamiltonian equations turn out to be 

{ 9 =Yg 

\ A = (ad(y))*A 

Let us now consider a linear vector held X, whose associated derivation is D = — ad(T’), and let 
us define the Hamiltonian 


H(X,g) = (Xg, Xg) = (A, TLg-lXg) = (X,Fg) 

According to Corollary [3] (in the Appendix) we have for any f£g 
B H 

— .Yg = (A, DY + ad (F g )Y) = ((D + ad^))* A, Y ). 

Finally the Hamiltonian equations for the Hamiltonian H = (A, F{g)) associated to the linear 
vector held X are: 


f 9 =* 9 

1 A =(D + ad (F g ))* A 


Remark: the inner case 

If the derivation D is inner, that is if D = —adX for some X £ 0 , then X has the form 
Xg = TLg.X - TRg.X and F(g) = X - Ad {g-^X. Consequently 

D + ad (F g ) = D + adX - ad(Ad(c/ _ 1 )A) = -ad(Ad( 5 ” 1 )A), 

and the second equation reduces to 

A = (-ad(Ad( 5 - 1 )A))* A. 


4.4 Hamiltonian equations of a simple ARS 

Consider an ARS defined as previously by (X, Yj .... ,Y n _ i), and consider the Hamiltonian 


/ n_1 \ i / n ~ 1 

n u (X,g,v,ui, . . ., u n —i) = ( X,vX X^UjYj) ) ~ + ^2 


u 


3 ' 


The associated equations are 


g = V X + YJI 1 UjYj 


A = (vD + ad ( vF(g ) + 1 u j Y j ) ) X 


In particular the equations of the normal extremals are obtained by application of the Pontryagyn 
Maximum Principle (see for instance [26], id, a, [25] i with v = 1. As usual the maximization of 
Ri w.r.t. v, u \,..., u n - 1 gives v = (A, X) and Uj = (A, Yj) for j = 1,..., n — 1. The maximized 
Hamiltonian is: 

1 1 n ~ 1 

H 1 {\g) = -{\X) 2 + 2 Y j {\Y j ) 2 . 

The case of the abnormal extremals is treated in Section [5] 


4.5 Semi-simple Lie groups 


In the semi-simple case all the derivations are inner and the remark of the end of Section 4.3 
applies. Moreover there exist on g an invariant scalar product. In order to avoid confusion it 
will be denoted by (. ,.} while the duality bracket will be denoted by (. , ,) d , and the invariance 
of the scalar product means 


VX, Y,Z eg (ad (X)Y, Z) a = - (Y, ad (X)Z) S . 

Thanks to this scalar product, we can identify g and g* by A' £g i —> A_y = ( X ,.) £ g*, so that 

(ad(X)*Ay, Z) d =(Ay,ad (X)Z) d 
= (Y,ad(X)Z) a 
= -(ad (X)Y,Z) a 

and ad(X)* can be identified with — ad(X). The previous equations become: 

f g = V X + UgYj 

\Z =[vM{g~ 1 )X-Yr 1 u j Y j ,Z} 


where Z £ g. 



5 Abnormal extremals 


Let (A t € [0, T], be an abnormal extremal. We know that 

(i) A(t) does not vanish; 

(ii) (X(t),Yi) = 0 for i = 1,... ,n - 1; 

(iii) (\(t) oTL g(t) -i,X g{t) } = (A (t),F g{ $) =0. 

These conditions imply that g{t ) belongs to the singular set Z for all t. Consequently the linear 
vector held X is a linear combination of Yf,..., Y n _\ along the curve g(t ) and we can assume 
the control v to vanish. This may modify the optimality of the associated control, but not the 
geometry of the curve, which is what we are looking at. 

Thanks to v = 0, the Hamiltonian equations reduce to: 

g = sr 1 «iYj 

* = (Ei _1 u M( Y jj) x 

As the one-form A (t) is not zero but vanishes on Y \,...,. Y n _\, that is on the left-invariant 
distribution A, it can be written as A (t) = p(t)u>, where p is a non-vanishing absolutely continuous 
function from [0, T] into M, and the ODE satisfied by A becomes 

p(t)u = p{t) 



This equality between one-forms is equivalent to the existence for almost every t € ([0, T]) of a 
real number c(t) such that 

(b) p(t) = c(t)p(t) 

Let us analyze the first equality. Let Y n G g such that uj(Y n ) = 1, and let Y £ A. Then 


n— 1 


ufdd(Yj) oj = c(t)uj 


3c G M s.t. ad {Y)*u = cu 

^ 3c G M, MX G g (w, ad(L') X) = c(u, X) 
•<=^- (oj, ad(y)li) = c(u), Yi) = 0 i = 1,. .., n — 1 
and {u, ad(Y)Y n ) = c{uo, Y n ) = c x 1 = c 
ad(T)A C A 


Let us set Y(t) = Ujjt)Yj. As g{t) is an abnormal curve, Y(t) belongs to A and satisfies 
ad(T(f))A C A for almost every t G [0, T]. In other words Y(t) satishes: 


Y(t) G Ap)jV'(A) a.e. t G [0,T], 
where jV(A) stands for the normalizer of A in g. Moreover 

c(t) = (uj, ad(y (t))Y n ) 

is measurable essentially bounded as soon as 1 1 —> Y ( t ) is. 


(8) 



Let us now look at the other condition, that is g(t) belongs to Z. As previously we denote 
by if the function from G to M defined by i/j(g) = (u g ,X g ) = ( 0 J,F g ). Recall that for all Y G 0 

Tg^.TLgY = (w, (D + ad (F g ))Y) . 

For Y G A P| AA(A) and g G Z we get 


Tgip.TLgY = (w,DY) . 


Indeed (uj, &d(F(g))Y) = —(cu,ad (Y)F(g)) = 0 according to F(g) G A and ([8]). Consequently a 
curve g(t) that verifies Q and g(0) G Z is contained in Z if and only if (u, Dg(t )) = 0 almost 
everywhere, that is if and only if g(t) G D~ l A almost everywhere. 


Conversely let g(t) be an absolutely continous curve defined on [0, T\ and such that g(t) = 
f t g{t) belongs to 


a = Af]M(A)f]D~ 1 A 


for almost every t. It is straightforward to check that it is the projection of an abnormal extremal 
(A (t),g(t)). The covector A (t) is equal to p(t)u, where p(t) is the solution of the linear equation 
p = c(t)p, with p(0) / 0 and c(t) = (w,ad (g(t))Y n ). 

We are now in a position to state the following theorem. 


Theorem 3 The vector subspace c\ of q defined by 

a = Af~]U(A)f]D~ 1 A 

is a subalgebra of g. 

It generates a connected Lie subgroup of G denoted by A, and for all g G Z the coset gA is 
included in the singular locus Z. 

The projections of the abnormal extremals are contained in Z and for g G Z the abnormal 
curves starting from g are all the absolutely continuous curves contained in the coset gA. 
Moreover the covecteur A (t) is up to a constant equal to p{t)u: where p{ 0) 0 and 

p(t) = (u>, ad(g(t))Y n ) p(t) 


where Y n is such that (oj, Y n ) = 1. 

Proof. We have only to prove that a is a subalgebra, and that gA is contained in Z as soon as 
g G Z. 

An element Y of g belongs to a if and only if 

(*) Y G A (ii) ad(T)A C A (in) DY G A. 

Let y, Z G a. By (ii) we get [Y, Z] = ad (Y)Z G A. 

Then ad(Y)ad(Z)A C ad(Y)A C A according to (ii), hence ad([Y, Z\)A C A. 

To finish D[Y,Z\ = [DY, Z] + [Y, DZ] G A since DY (resp. DZ) belongs to A and ad (Z) 
(ad(Y)) sends A into A. 

Consequently [Y, Z\ belongs to a which is a subalgebra of g. 

To prove the second point let g G Z. Any absolutely continuous curve starting from g 
and contained in gA, in particular any curve of the form gexp(tY) with Y G a, satisfies the 
conditions discussed before the statement of the theorem. Consequently it is an abnormal curve, 
it is contained in Z. 

□ 


Remarks 


1. The dimension of a is at most n — 2. Indeed a is included in A and its dimension cannot 
exceed n — 1. Suppose it is equal to n — 1, then A would be a subalgebra (because it would 
be included in J\T (A)), and would be moreover included in D~ 1 A, hence invariant for D. 
But if A is a D-invariant subalgebra then the rank condition does not hold. Consequently 
dim(a) < n — 2 . 

On the other hand this dimension is equal to n — 2 as soon as A is a subalgebra since in that 
case a = A P| D~ 1 A, and A cannot be included in D~ 1 A, which is therefore a codimension 
one subspace of g distinct from A. 

2. If F(g) never belongs to the subalgebra a (at no g E Z or at least at no g' E gA) then the 
controls associated to the abnormal extremals are uniquely defined and cannot involve the 
linear vector field. Consequently the abnormal curves are the left translation by g of the 
geodesics of A for the left-invariant metric induced by the one of the ARS. 


6 Almost-Riemannian and sub-Riemannian structures 

In this section we show that a simple ARS on a n-dimensional Lie group is related to a codimen¬ 
sion one regular sub-Riemannian structure on a (n + 1 (-dimensional Lie group. In other words 
the desingularization of a simple ARS is global and the desingularized structure is a classical 
sub-Riemannian one. 


6.1 The lift of a linear vector field 

Let G be a connected Lie group, g its Lie algebra identified with the set of left-invariant vector 
fields, and X a linear vector field on G. The flow (<^t)teM of A is a one-parameter group of 
automorphisms of G, hence a Lie group morphism from M into G, so that we can consider the 
semi-direct product G = G x ^ M defined by the law 

(5l,n)-(52,T 2 ) = (</V 2 ( 51 ) 52 , Tl +T 2 ) 

The group G can be identified with the closed and normal subgroup G x {0} of G, and M with 
the closed subgroup {e} x M of G. This last is normal if and only if X = 0 (see [ 8 ] for that 
construction). 

Let g stand for the Lie algebra of G. It is equal to some semi-direct product g x K ; in view 
of its characterization the action of the Lie algebra M on g should be identified. Let X denote 

0 _ 

the left-invariant vector held which is equal to (0,1) = —— at the identity e = (e, 0). This held 

~ de¬ 

generates the Lie algebra M, so that exp (tX) = (e,t). As it is left-invariant, we get also: 

V( 5 , t) e G ( 5 , r) exp (tX) = ( g , r).(e, t) = {y t {g),t + r), 


and 

T(5,t) 

Thus we get the equality X = X + 
D: 


= 4 . (Vt{g),t + T) = {X(g), 1). 

at \t=o 
d 

——, and denoting as usual the derivation associated to X by 
OT 


VF E g [X, Y] = [X, Y] = -DY. 
Finally the Lie algebra g is the semi-direct product g x_£> M: 


[*1 + nX, Y 2 + t 2 X] = [Li, Y 2 \ - t\DY 2 + t 2 DY x . 


6.2 Projection 

The group G is diffeomorphic to the homogeneous space G/M, the set of right cosets of M, but 
this last is not a quotient group in general because the subgroup M of G is not normal as soon 
as X ^ 0. On the other hand the left-invariant vector fields of G can be projected to G/M (this 
is always true for left-invariant vector fields and right cosets). 

More accurately we get 

( e ! r ) exp {t(Y, 0)) = (e,R)(gexp(tY),r) = (gexp(tY), M) 
(e,M)(ff,r)exp(tT’) = (e,R)(<p t (g), r + t) = (<p t (g),R) 

so that we can identify the projection onto G/M of the vector field (Y, 0) (resp. X) with the 
left-invariant vector field Y (resp. the vector field X) of G. 

6.3 Almost-Riemannian and sub-Riemannian structures 

Let us now consider an almost-Riemannian structure defined on G by 
{X, Yi,, Y n _ i}, assumed to satisfy the rank condition, which implies 

CA{X,Y 1 ,...,Y n _ 1 }=RX®g. (9) 

On G, we can define the left-invariant distribution 

A = Span{A,y 1 ,...,y n _ 1 } 

This distribution is left-invariant, of codimension 1, and on account of ([9| we get 

CA{X, Ti,..., y„_i } = RX ©0 = 0 

which shows that A satisfies the rank condition. In order to define a corank one classical sub- 
Riemannian structure on G it remains to declare the vector fields X, Y \,..., Y n -\ to be orthonor¬ 
mal. 

Notice that the dynamics of the almost-Riemannian structure being g = vX + u\Y\ + ■ ■ ■ + 
v n -{Y n _i the one of the sub-Riemannian structure writes in the coordinates ( g , r) G G xi M: 

f 9 =vX + UlYi H-h v n -\Y n -i 

\ t = v 


6.4 Optimality 


Let u stand for the general control (v, u \,..., u n - 1 ), let go and g\ be two points of G, and let 
ti be a minimal control that steers go to g\ on the time interval [0, T]. By a minimal control is 
meant a control such that the associated curve 7 (f) minimizes the length between go and gi, and 
such that 

71—1 


“(Oil 


N 


v 2 tf) + = ii'vwn 

%— 1 


a.e. t G [ 0 ,T], 


On G the control u steers 


(.90i 7"o) to (gi, t 1 ) for any tq and T\ that satisfy 


n - r 0 


v(t) dt 


( 10 ) 


and it is minimal. 




Indeed let tGK, let u be a minimizing control that steers (go,0) to (gi,r) in time T, and 
let 7 be the associated trajectory. Its length is: 

l( l) = f \\u(t)\\dt. 

Jo 

Let 71 be the projection of 7 on G. As it steers go to gi, we get: 

[ ||u(t)|| dt = 1 ( 7 ) < Z( 7 i) < l{j) = [ \\u(t)\\dt. 

Jo Jo 

In conclusion the extremal trajectories on G can be lifted to extremal trajectories on G: the 
optimal control u belongs to the set of optimal controls that steers (go,0) to (<?i,t) for some r, 
but this last is not known "a priori". 


7 Examples 


7.1 An example on the 2D affine group 


With the notations of Section 3.5.2 we consider the simple ARS defined by the left-invariant 
vector field X and the linear one X = Y +I*Y, that is: 



0 x — l\ 

0 0 )■ 


The pair {X,X} is considered as an orthonormal frame, which defines the almost-Riemannian 
metric. Notice that 


(i) these two vector fields are linearly independant on the open and dense subset { 1 ^ 1 }; 

(i) the rank condition is satisfied since [A, A] = Y. 

Consequently the ARS is well defined and its singular locus is the line Z = {x = 1}. 

The dynamic of that ARS is described by g = vX + uX, where v, u £ R, or in coordinates: 

(e> 

The associated Hamiltonian is 

H = upx + vq(x — 1) — - 
where (p. q) stands for the covecteur A. 


ux 

v(x — 1 ) 


2 


,, 2 \ 


7.1.1 Abnormal extremals 

We know by Theorem [3] that the abnormals are contained in some cosets of a subgroup of 
dimension at most n — 2 = 2 — 2 = 0. Let us verify that on A// + (2) the abnormals are fixed 
points. 

If v = 0 we get 

Hq = upx + vq(x — 1) = u (A, X) + v (A, X) 

so that Ho can have a maximum w.r.t. ( u , v) if and only if (A, X) = (A, X) = 0 (here (.,.) stands 
for the duality bracket). As A / 0 the fields X and X should be linearly dependant, which 
happens only in the singular locus. But then y = 0 and the abnormal extremals are fixed points. 



7.1.2 Normal extremals 


As usual, v can be normalized to v = 1. The computation of the partial derivatives of T~L w.r.t. 
u and v shows that T~L reaches its maximum for u = px and v = q(x — 1 ), so that this maximum 
is equal to 

H = -p x + -q (x - 1) , 

and the Hamiltonian equations are: 

J x= px 2 J p= —p 2 x — q 2 (x — 1 ) 

\ y= q(x-l) 2 \ <7 = 0 

In the particular case where q = 0 we get p 2 x 2 = 2 H, and setting c = \/2 H we obtain 
x = ±cx and y = 0 . 

Consequently the lines parametrized by x(t) = xoe± ct and y{t) = yo are geodesics. 


7.1.3 Poincare coordinates 


We are interested in the geodesics starting from the singular locus, that is with initial conditions 

Oo = 1 ) yo)- 

We can restrict ourselves to the geodesics parametrized by arclength, that is to H = —. As 

H = ~pq at t = 0 for geodesics starting from (1, yo), this amounts to take po = e = ±1. Moreover 
we can also assume q = qo > 0 because of the symmetry w.r.t. the x axis. 

To solve the Hamiltonian equations two changes of variables are done. Since H is constant 
we can first set: 

ecos(a) = px 
esin(a) = q(x — 1 ) 

Derivating the first of these formulas we get 

—ea sin(a) = px + px = p 2 x 2 — p 2 x 2 — q 2 x(x — 1 ) = —qxe sin(a), 

so that 

a = qx = esin(a) + q. 

For an initial point in the singular locus, that is xq = 1, and according to q 7 ^ 0, we get ao = 0 
(or kir). The equation to be solved is consequently 


a = po sin(a) + q with cc(0) = 0 and q > 0. 


( 11 ) 


,a\ 


In order to solve (11) we do the change of variable r = tan( —), and the equation becomes: 


r = er + -q(l + r 2 ) r(0) = 0. 


(12) 


Before solving (12), we can notice that a straightforward computation gives: 

e 2 t 

X = 1 + 


ql + T 2 


e 2 t 2 

y = qt + - ——~ 2 arctan(r) + y 0 
q 1 + t 1 







The separation of the variables in (12) leads to 


2 I" dr 

q J t 2 + 2 rr + 1 


where r = -. To compute these primitives we should distinguish three cases according to the 

q 

sign of the discriminant 4 (r 2 — 1) of r 2 + 2rr + 1. 

Case 1 (0 < q < 1) 

We get 


9 (ex P (^? t) - l) 


r(t) = 


We are interested in the first return to the singular locus {x =1}. It happens when r(t) 


+ y/l - q 2 + ( v 7 ! - q 2 - e) exp( v / l - q 2 t ) 


for some to, or as a limit if r(t) 


±oo. 


For 0 < q < 1 the function r(f) does not vanish for t > 0 but tends to +oo when t 




tends to 


: In 


: + \/l 

- 


which is positive only for e = 1. Consequently when e 


\/l — 


: In 


1 + 

1- 


, then x(t) tends to 1, and 


vanishes 

tends to 
1 and t 


y(t) 


q 


In 


1 + y/l~g 2 \ 

i-\/i — q 2 J 


2 

H-7T + 2/0- 

q 


Notice that the geodesics do not go back to the y-axis in positive time for e = — 1, that is for 
x < 1, as illustrated on the following picture: 


Geodesics for q=0.5 



Case 2 (q = 1) 
We get r(t) 


- hence 

2 -et 


x(t) = 1 + 


y(t) = t + 


2 et - t 2 
2-2 et + t 2 

et 2 

2-2 et + t 2 


2 arctan( — ) + y 0 

£ 6C 




































The trajectory belongs to the singular locus, that is x(t) = 1, for t = 0 and t = 2e. 

When 1 1 —> 2e, y(t) tends to e(4 — ir), that is the locus of the first return to the y-axis. Notice 
that as well as in the first case the geodesics do not go back to the y-axis in positive time for 
e = —1, that is for x < 1. 



0 0.5 1 1.5 2 


Case 3 (q > 1) 
We get 


where 6 = arctan 


r(f) = —^ tan ^-£^/ q 2 — 1 + ed^j — 


1 


yV _l / 

As in the hrst case r(£) does not vanish, hence x(t) cannot be equal to 1, for t > 0, but 
r(f) i —y Too when t i-t — t Consequently when t tends to that limit we get: 


x{t) = 1 T 


e 2 r 


q 1 T t 2 


e 2 r 


y(£) = qt T - —— 9 - 2 arctan(r) T yo '—> q 
q 1 T t z 


it - 2ed 2e 


Vq 2 - i ^ 


H-vr + yo 


The behaviour is here different from the two previous cases because the geodesics return to the 
axis {x = 1} in both cases e = ±1. 


















Geodesics for q=2 



To finish the three behaviours are compared in a unique picture and the unit sphere is drawn 

Geodesics Unit Sphere 




7.2 An example on the Heisenberg group 


The equations of the following example are only stated, not solved. Actually the purpose is here 
to illustrate some of the theoretical results, in particular the desingularization process. 

With the notations of Section |3.5.3 we consider on the Heisenberg group G the ARS defined 
by the set of vector fields {A, A, Z} where 


v <9 1 9 d 

X = X By + 2 X Bi 


is associated to the derivation D 


/0 0 0 \ 
1 0 0 
\0 0 0 / 


Here A = Span (A, Z } is an ideal of g, and according to Proposition [I] the singular locus 
is a subgroup. An easy computation shows that it is equal to the set of fixed points of A: 
Z = {x = 0}. On the other hand the rank condition is satisfied since [A, A] = DX = Y , and 
the ARS is well defined. 

The dynamic of that ARS is here g = vX + u\Y\ + U 2 Y 2 , where v, U\,U 2 £ M, or in coordinates 


(£) 


X = U\ 
y = vx 

Z = U2 + \VX 2 


and the associated Hamiltonian is 

% = (A, vX + U 1 Y 1 + U 2 Y 2 ) - \u{v 2 +u\ + u |) 

= pu\ + qvx + r{u 2 + \vx 2 ) — \v{v 2 + u 2 + n|), 


v > 0. 





























where A = (p, q, r ). 


7.2.1 Abnormal extremals 

Let us compute the abnormals extremals. For r = 0we get 

%o = U \P + u 2 r + v(qx + ^rx 2 ) 

so that T~Lq presents a maximum if and only if p = r = qx = 0. But p = r = 0 imply q ^ 0, hence 
x = 0. Moreover q is a constant because — o = 0, and choosing g=lwe finally obtain x = 0 
(that is the abnormals are contained in the singular locus) together with the equations 

f y =o 

\ z = u 2 

In other words the vertical lines through a point (0,yojO) are abnormal extremals. 

On the other hand the normalizer of A = Span{X, Z} being g and D _1 A being equal to 
Spanjy, Z}, we know by Theorem [3] that the abnormals extremals are the curves contained in 
Z = {x = 0} and parallel to a = Span{Y}, that is the computed result. 


7.2.2 Normal extremals 


As previously v is normalized to v = 1. The computation of the partial derivatives of Tii w.r.t. 
v, u\ and u- 2 shows that TL\ reaches its maximum for u\ = p, u 2 = r , and v = qx + so that 
this maximum is equal to 



p 2 + r 2 + (qx + -xx 2 ) 2 


and the Hamiltonian equations are: 


x = p 

y = qx 2 + \rx 3 
z = r + ^qx 3 + |rx 4 


p = — (q + rx)(qx + |rx 2 ) 

q =o 

r =0 


One can see at once that the abnormal extremals are normal, they are obtained for x = p = 0. 
On the other hand the case q = r = 0, p / 0 gives x = p and p = y = i = 0, hence all the 
lines parallel to the x-axis. 


7.2.3 Poincare coordinates 

Apart from these simple cases, we are lead to isolate the two following equations: 


( x = p 

\ p = — (q + rx)(qx + \rx 2 ) 


1 


1 


2\2 


We know that H = - + r 2 + qx + (qx + -rx 2 ) 

as r is also constant, we can set 


is constant along the extremals and, 


and then 


2 , / I 1 2\2 n rj2 2 2 

p + (qx + -rx ) = 2 H — r = c , 

J c cos a = p 
\ csina= qx + \rx 2 


By a first derivation we get 


J — casino; = p =—(q + rx)(qx + ^rx 2 ) =—(q + rx)c sin a 
| ca cos a = (q + rx)x = (q + rx)p = (q + rx)c cos a 

so that a = q + rx. A second derivation gives: 

a = rx = rp = cr cos ct, 


that is the pendulum equation. 

For an initial point in the singular locus, that is xo = 0, we get ceo = 0 (or fc7r), and po = ±c. 
The equation to be solved is consequently 

a = p 0 r cos(a) with a(0) = 0 . 


7.2.4 Desingularization 


The Lie algebra generated by A, Y\. Y 2 is equal to Sp{X, X, Y, Z} = g©RA. This algebra is the 
Engel one since [A, A] = Y. [X. Y] = Z, and the other brackets vanish. According to Section [6] 
it it isomorphic to the Lie algebra of the semidirect product of G by M, that is g © RA where 
~ d ~ 

A = A + ——. The lift of the almost-Riemannian distribution to G = G x R is the dimension 
aw ~ 

3 classical sub-Riemannian distribution A = Span{A, Yi, I 2 }, the dynamics of which writes in 
natural coordinates: 

x = U\ 
y = vx 

Z = U2 + \vX 2 
W = V 

The associated Hamiltonian is 


% =/ H + sv 


v > 0. 


where A = (p, q, r, s) £ g*. 

Abnormal extremals. An easy computation shows first that v = 0 implies p = r = s + qr = 0, 
and then that the abnormal trajectories are the lines parallel to Oz (after a distinction between 
the cases s = 0 and s / 0). 

Normal extremals. Setting v = 0 we obtain at once the maximized Hamiltonian 



+ r 2 + (s + qx + \rx 2 ) 2 ) 


and the Hamiltonian equations: 
x = p 

y = sx + qx 2 + T^rx 3 J p = ~(q + rx)(s + qx + \rx 2 ) 

z = r + \sx 2 + ^qx 3 + |rx 4 ( q = r = s = 0 

w = s + qx + ^ rx 2 

The almost-Riemannian dynamics is the part of that last obtained for s = 0. 




8 Extension of the model and Equivalence 


The model analyzed in the previous sections can be generalized in two ways. Firstly on a 
connected ra-dimensional Lie group G, instead of considering n — 1 left-invariant vector fields 
and one linear, we can consider a set {F\,..., F n } of n affine vector fields. By affine vector held 
is meant an element of the normalizer of the Lie algebra g of G in the Lie algebra of analytic 
vector fields on G. The affine vector fields are sums of right-invariant, left-invariant and linear 
vector fields (see [22]). so that they are either left-invariant or the sum of a linear vector held 
and a right-invariant one. 

One of the reasons for considering affine vector helds instead of linear ones is that it may 
happen that the zero locus of several non left-invariant vector helds contain no common point, 
which would be the case, the common point being the identity e, if all these helds were linear. 

Definition 2 A general almost-Riemannian structure on a connected n-dimensional Lie 
group G is defined by a set (Fj,..., F n } of n affine vector fields that satisfy the two following 
properties: 

(i) The rank of the distribution defined by Span{F ±,..., F n } is equal to n on a proper open 
subset U of G. 

(ii) This distribution satisfies the rank condition. 

The almost-Riemannian metric is defined by declaring the set {Fj,..., F n } of vector fields to be 
an orthonormal frame. 

The second generalization consists in considering the same framework on homogeneous spaces. 
Let H be a closed subgroup of G, and G/H the manifold of right cosets of H. It is a well known 
fact that the left-invariant vector helds can be projected to the quotient group. The set of their 
projections will be denoted by LLg; it is a Lie algebra. On the other hand an affine vector held 
on G/H is defined in a natural way as the projection, whenever it exists, of an affine vector held 
of G. However it is shown in |22] that the affine vector helds of G/H are exactly the elements of 
the normalizer of n*g in the Lie algebra of analytic vector helds on G/H. Consequently we can 
copy the previous definition: 

Definition 3 A general almost-Riemannian structure on a connected n-dimensional ho¬ 
mogeneous space G/H is defined by a set {Fj,..., F n } of n affine vector fields that satisfy the 
two following properties: 

(i) The rank of the distribution defined by Span{F\,... ,F n } is equal to n on a proper open 
subset U of G/H. 

(ii) This distribution satisfies the rank condition. 

The almost-Riemannian metric is defined by declaring the set {Fj,..., F n } of vector fields to be 
an orthonormal frame. 

Remark. In both cases all the vector helds are analytic and the condition that the rank of the 
distribution dehned by Span{Fj,..., F n } be equal to n on an proper open subset U implies that 
U is dense in G (or G/H). 

Consider now a connected and n-dimensional manifold M endowed with an almost-Riemannian 
structure dehned by a set {/i,..., f n } of vector helds that verihes 

(i) The rank of Span{/i,..., f n } is equal to n on an open and dense subset U C M. 


(ii) The set {/i,..., f n } satisfies the rank condition. 

As previously the almost-Riemannian metric is defined by considering the vector fields {f\,..., f n } 
as an orthonormal frame. 

If we want this structure to be equivalent to a general almost-Riemannian structure on a Lie 
group or a homogeneous space, we should add the two following necessary conditions: 

(iii) The vector fields f n generate a finite dimensional Lie algebra. 

(iv) The vector fields fi, ■.., f n are complete. 

These conditions turn out to be sufficient. 

Theorem 4 Under the conditions (i) to (iv) the manifold M is diffeomorphic to a homogeneous 
space G/H and the vector fields f\,..., f n are related by this diffeomorphism to invariant or 
affine vector fields on G/H so that the almost-Riemannian structure on M is equivalent to a 
general almost-Riemannian structure on G/H. 

If condition (iv) is removed, then the statement holds locally. 

Proof. Let £ be the (finite dimensional) Lie algebra generated by {/i,..., f n } and for each subset 
/ of {1,2,... ,n} let us denote by £(/) the ideal generated in £ by {/*; i 6 I}. Let us choose 
/o as one of the subsets of { 1 , 2 ,..., n} that satisfy the rank condition and whose cardinal is 
minimal among the subsets of {1, 2,..., n} that satisfy that condition. This is possible since on 
the one hand £ = £({ 1 , 2 ,..., n}) satisfies the rank condition, and on the other hand the choice 
is done in a finite set. 

We set £o = £(/o)- Let G be the simply connected Lie group whose Lie algebra g is isomorphic 
to £o- According to the Equivalence Theorem of |22j there exist a closed subgroup H of G and 
a diffeomorphism 4> from G/H onto M such that 

£o = ^(LRg), 

where 11*0 stands for the Lie algebra of invariant vector fields of G/H. 

We can assume without loss of generality that Iq = {1,... ,k} with k < n. As £o is an 
ideal of £, the fields fk+i, ■ ■ ■, fn belong to the normalizer of Cq. Their images Fi = 
i = k + 1 ,..., n belong to the normalizer of 11 * 0 , hence are affine. 

If we do not assume the ff s to be complete then the Equivalence Theorem applies locally. 
This is not explicitly stated in [22], but is clear from the proof of the Equivalence Theorem 
(Theorem 5.1 of [22]). 

□ 


Corollary 2 Under the previous conditions, if moreover M is simply connected and dim(£o) = 
dim(M), then the ARS on M is diffeomorphic to a general ARS on a Lie group. 

If M is not simply connected, or if the vector fields fi are not complete then the same statement 
holds locally. 

Proof. With the notations of the previous proof we get: if dim(£o) = dim(M) then dim(M) = 
dim(G'), the subgroup H is discrete, hence reduced to the identity because of the simply con¬ 
nectedness assumption on M. 

□ 

Remark (Due to Jean-Paul-Gauthier) 

In the statement of Theorem [4] the ARS on the manifold M is globally defined by n orthonor¬ 
mal vector fields that generate a finite dimensional Lie algebra. The problem to know if a general 
ARS can be defined, at least locally, by a finite set of orthonormal vector fields that generate a 
finite dimensional Lie algebra is open. 


9 Appendix 


Let X be a linear vector field on a connected Lie group G. We denote by F the mapping from 
G into 0 defined by 

F{g) = TL g -iXg (13) 


Proposition 4 1. For all Y e 0 


dt \t=o 


F(exp(tY)) = DY 


(14) 


2. For all g 6 G and for all Y £ 0 

F(g exptY) = F(exptL') + e~ tad ^ F(g) (15) 

and 

T g F = (D + ad(F g )) o TL g - 1 (16) 

Proof. 

1. Denoting by (<pt)teR the flow °f X and according to Formula ([3| (see Section [2]) we have: 

dt|t=o iJ ’ (eXp( * y)) = Jt\t=o TL{eMtY)) ~ lXexp{tY) 

= i y- exp (—tY)ip s (exp (tY)) 
dt\t=o ds\ s =o 

= -f- -f- exp(—ty) exp(te si:, y) 

ds|s=odf|t=o v ; 

= (-y + e sl? y) = py 

ds|s=o 


2. Let y e 0 . According to Formula [2] of Section [2] we get: 

F(gexp(ty)) = TL(g exp(ty ))- 1 Xg exp (tY) 

= T L exp (_ tY )T Lg_l(TL g X exp ( tY ) + TR eMtY) Xg) 

= F(exp(ty)) + Ad(exp(— tY))F(g) 

= F(exp(tY)) + e- tad{Y) F(g ), 


that is Formula (15). Derivating that formula at t = 0 we obtain at once Formula (16): 


T,F.Y g = TgF.TLg.Y = | |( _/(9exp(tr)) 

= s |1 . 0 ( F ( exp < (r » +e "“ d(i ' lF < ! ' ) ) 

= DY - ad (Y)F(g) = DY + ad (F(g))Y 


The next very useful corollary is immediate from Formula (16). 


□ 


Corollary 3 Let uj be a left-invariant one-form, on G, and if the real valued function defined on 
G by ip(g) = ( uj,F(g )). Then 


vy € 0 Tgif.Yg = (uj, DY + ad(Fg)Y). 
This can also be written T g ijj °TL y = (DY + ad(F g ))* uj. 





To finish we compute F as a power series. 

Proposition 5 For all g E G, Y E g, and t E M: 

/jk 

(i) Vk > 0 ^F{gexptY) = (-l^a/"^^ + (-l) fc ad fc (F)F( 5 exptF) 

+°° yfc 

(**; T(exptY) = E (-l)--ad-(P)W 
fc=l 


Proof. 


(i) According to Formula (16) of the previous proposition we get 


d 

dt 


F(gexp(tY )) = T gexp(tY) F.Y gexp(t Y) = DY - ad(Y)F(g exp(fY)), 


that is the expected formula for k = 1. By induction we have also: 


(Jk-\-l 

—F(gexptY) = (-l) fc ad k (Y)(DY - ad(y)F( 5 expty)) 

= (-l) fc ad fc (y)Ziy + (—l) fc+1 ad fc+1 (y)F(g exp tY). 


(ii) We have y(exp(0y)) = F{e) = 0 and for k > 1 the formulas of item (i) applied to g = e 
and t = 0 gives 

|^ =( F(exp(fy)) = (-l) fe - 1 ad fc “ 1 (y)T)y, 

+°° t k 

so that by analycity A (exp tY) = y^(—l) fc ~ 1 — ad fc ~ 1 (y)Ziy. 

k= 1 

□ 
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